Peak Effect versus Skating in High Temperature Nanofriction 
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The physics of sliding nanofriction at high temperature near the substrate melting point Tm is so 
far unexplored. We conducted simulations of hard tips sliding on a prototype non-melting surface, 
NaCl(lOO), revealing in this regime two distinct and opposite phenomena for plowing and for grazing 
friction. We found a frictional drop close to Tm for deep plowing and wear, but on the contrary 
a frictional rise for grazing, wearless sliding. For both phenomena we obtain a fresh microscopic 
understanding, relating the former to "skating" through a local liquid cloud, the latter to linear 
response properties of the free substrate surface. It is argued that both phenomena should be 
pursued experimentally, and much more general than the specific NaCl surface case. Most metals 
in particular possessing one or more close packed non-melting surface, such as Pb, Al or Au(lll), 
that should behave quite similarly. 



Thanks to the advent of increasingly sophisticated tip- 
based tools, the physics of nanofriction is now a field of 
growing importance. There is in particular a need for a 
deeper understanding of the physics of friction in extreme 
regimes including very high temperatures (HT) where the 
sliding materials approach their melting points. In ev- 
eryday life, high frictional temperatures are routinely at- 
tained in a variety of situations including machining, mo- 
tors, braking, etc While rather complex phenomena 
are generally involved in many such practical situations, 
at the fundamental level even the very basic microscopic 
friction processes taking place close to the substrate melt- 
ing point are not sufficiently explored, characterized and 
understood. For bare metals on metals, for example, a 
high speed drop of wear-related friction is generally re- 
ported [IJ , and reasonably attributed to softening caused 
by high ("flash") temperatures ("if one could move fast 
enough, one could ski on copper mountains. . . " [2J), but 
not yet pursued with surface science standards. Con- 
versely, dry grazing friction between hard materials has 
been known to rise at HT, ostensibly due to the onset of 
some diffusion [2] ; again however that has not been fur- 
ther followed up microscopically. At the present stage, 
it is just not clear when, how and why closeness to the 
substrate melting point will imply a drop or a rise of 
frictional forces under controlled contact circumstances. 
Nanofriction presently offers a fresh opportunity to ad- 
dress the microscopic physical facts behind these phe- 
nomena which are receiving increasing attention. 0] Re- 
cent simulation studies have addressed metallic friction 
between a hard indenter and a soft surface, at high veloc- 
ities and at tip-driven local temperatures near and above 
the bulk melting point [4(. Also the simulated motion 
of a rigid tip across a Si surface at room temperature 
showed an interesting tip and temperature driven local 
phase transition leading to a reduced friction [5| . 

Considering that most AFM nanotips slide normally 
at very low speed, and can hardly perturb the massive 
substrate, our attention is focused on friction under con- 




FIG. 1: Left panel: rigid diamond tip constructed for plowing 
friction. The tip length is ~ 26 A, the diameter is about 
13 A. Right panel: rigid diamond tip constructed for wearless 
friction, (a) side view; (b) top view. The tip length is ~ 7 A, 
the diameter is about 13 A. 



ditions of essential closeness to thermal equilibrium, how- 
ever at temperatures arbitrarily near or even marginally 
above the substrate Tm- We arrive at a first charac- 
terization of possible scenarios by means of simulation 
and theory. Our starting point are molecular dynamics 
(MD) simulations of sliding nanotips on the hard crys- 
talline surface of a solid close to its bulk melting point 
Tm- We actually recover both the high temperature fric- 
tion drop for deep plowing and wear, and the frictional 
rise for weak load wearless sliding; and for both regimes 
we obtain a fresh microscopic understanding. 

For a meaningful case study we require first of all a 
careful choice of surfaces and tips. In fact, close to Tm 
most solid surfaces wet themselves - in full equilibrium 
conditions and without any perturbing tips - with a mi- 
croscopically thin film of melt (surface melting) Q . The 
liquid film will generally "jump to contact" with the tip 
long before closest approach - to some extent ruining 
the experiment, as indeed observed Q and predicted [s| 
long ago. This unwanted complication could be avoided 
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for example by resorting to tips that are not wetted by 
the substrate. Alternatively one could, as we do in the 
present work, choose a hard, non-melting crystalline sub- 
strate facet exhibiting no spontaneous surface melting at 
all - a surface that will spontaneously remain solid and 
crystalline up to Tm @]. One added advantage of this 
choice is the possibility to distinguish between very dif- 
ferent frictional regimes, such as deep plowing (involving 
plastic flow and heavy wear) and superficial (wearless) 
grazing friction, a distinction otherwise made impossible 
by surface melting. We did not consider inhomogeneous 
surfaces, simply kept from melting by some kind of pro- 
tective coating, for they will be very difficult to control. 
Close-packed metal surfaces, such as Pb, Al or Au(lll) 
are known to be nonmelting [f| and would serve the pur- 
pose. Eventually, the very stable alkali halidc NaCl(lOO) 
surface proved close to an ideal choice because of its 
extreme and recently understood nonmelting habit @. 
One added bonus of alkali halide surfaces is that they 
are already choice substrates for abundant experimental 
nanofriction work at room temperature [Toj ] . 

The high temperature behavior of NaCl was simulated 
by the classic Born-Mayer-Huggins-Fumi-Tosi (BMHFT) 
two body potential [111]. Bulk molecular dynamics (MD) 
simulations [tj HH, as well as thermodynamic integra- 
tion yield for that model a bulk melting temper- 
ature Tm =1066 K (experimental value 1074 K). The 
NaCl(lOO) surface was studied using periodically re- 
peated slabs consisting of 12 24 planes of 100 NaCl 
units each separated by 100 -j- 120 A of vacuum. Compu- 
tational details were similar to those of Ref. [T2] The per- 
fect, tip free NaCl(100) was found to remain indefinitely 
solid and free of diffusion at all temperatures up to Tm 
and, in a metastable state, even above Tvr up to a "sur- 
face spinodal temperature" Tg ~ 1200 K, for our longest 
simulation times ~ 200 ps. This qualifies NaCl(100) as 
a nonmelting surface, in agreement with its very poor 
wetting by molten NaCl [(f 
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To simulate nanofriction we modeled a rigid diamond 
tip apex, for which deformations are negligibly small. For 
indentation and deep wear (plowing) simulations we used 
the sharp ~ 400 atom conical tip (Fig.QJ,). For wearless 
(grazing) frictional simulations we used instead a blunt 
flat tip whose contact plane is composed of ~ 200 atoms 
(Fig.[l]b) forming a diamond (111) plane of ~ 13 A diam- 
eter. We slid our tips at constant velocities and purposely 
avoided including any cantilevers with their mechanical 
deformations, and related stick-slip phenomena. Thanks 
to the relatively minor perturbation represented by the 
sliding tip, even in the plowing case no additional heat 
sink other than temperature control through standard 
velocity rescaling was necessary to ensure efficient dissi- 
pation of the Joule heat. 

We address first of all plowing friction, with the aim 
of understanding wear forces and their dependence on 
temperature. For that, we indented the surface with the 
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FIG. 2: Instantaneous forces in a plowing friction simulation 
of NaCl(100) under constant penetration of 6 A, speed of 
50 m/s, and T = 300 K 
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FIG. 3: Averaged plowing frictional force (F x ) as a function 
of temperature. Note the drop at T 3S ~ T m — 150 K. The 
frictional stress a — {F x )/A can be obtained using an active 
tip plowing area of ~6 A 2 . 



sharp tip and slid it at constant velocity v (generally be- 
tween 25 and 50 m/s) parallel to the NaCl(100). Because 
under constant load at high temperatures the tip pene- 
tration will generally be unstable with time, we worked 
in a constant penetration mode, the tip wedge entering 
the NaCl(100) surface by a fixed depth d = 6 A. The 
forces acting on the tip were collected during the sim- 
ulated plowing along the x direction - parallel to the 
surface - across the 60 x 60 A 2 wide surface. 

The vertical and in plane forces F z (t) and F x (t) are 
shown in Fig. [2] The large noise reflects the discontinuous 
nature of wear, including plastic flow and the formation 
and displacement of debris. As a function of temperature, 
the mean frictional force (F x (t)), shown in Fig. [3] clearly 
shows two regimes. 

Up to T = - 900 K ~ 0.85 T M , the plowing 
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frictional stress is approximately constant, correspond- 
ing to a friction coefficient \i = (F x /F z ) ~ f. This 
large value is naturally attributed to plastic deformation, 
breaking and reforming of atomic bonds, and the digging 
of a trough. For macroscopic indenters onto larger size 
substrates, plowing would in addition generally cause dis- 
locations to appear. Our nanosized, poorly indented tip 
is in fact unable to generate dislocations and only leads 
to local plastic deformation and eventually local melting 
very near Tm. The typical magnitude of F x (t) is around 
1 nN, which given the small tip active area is equivalent 
to a stress of about 6 GPa, a value close to the estimated 
yield stress of 8.3 GPa of NaCl against (100) shear near 

An approximate frictional energy balance at 300 K 
shows that about 20% is spent to plastically deform and 
scratch the material, and the remaining 80% is dissipated 
as heat, similar to typical figures given in the litera- 
ture [l5|. The work-induced local melting of NaCl(100) 
near the tip occurs at essentially constant temperature, 
the mechanical energy transfered into heat of melting, 
plus phonons that are effectively absorbed by the temper- 
ature control through velocity rescaling. Above T so ft ~ 
Tm— 150 K, the frictional force drops and remains low 
across the Tm and above up to T ss ~ 1200 K where the 
surface spontaneously melts. The drop of friction above 
T so ft resembles the experimental one reported for metals 
on metals, and is reminiscent of that observed for hard 
sliders on ice. In ice [TBI, the accepted view is that the 
low skating friction is due to a thin sliding-induced lu- 
bricating liquid layer. In our simulation, we find that 
indeed the frictional drop is accompanied by a thin local 
cloud of displaced Na and CI ions which surround the tip 
with liquid- like mobilities (see Fig. 3]). The local melting 
is motion-induced, and not due e.g., to uniaxial stress- 
induced lowering of the melting point. 0, E3] Coated 
by this nearly liquid shroud (Fig. [4|) the tip effectively 
skates. The ease with which the tip wades through the 
solid is remarkable. The scratch produced by the moving 
tip, permanent at low temperatures, heals rapidly away 
at higher temperatures, (see movie in the supplementary 
material) with a characteristic time which is in the ps 
range at Tvi- The fast resolidification behind the tip re- 
stores the equilibrium solid nature of the NaCl surface 
at Tm d, [HJ j and the resulting latent heat is dissipated 
away as phonons. 

Totally different from plowing friction are the results, 
shown in Fig. for wear less friction of a large flat tip, 
grazing the NaCl(100) surface under zero load. Here 
the tip contacts the substrate by mere adhesion, the 
adhesive energy V4d — 0.7 eV equivalent to an effec- 
tive load F z ~ — Vad/h where h is of order of half the 
interplanar spacing (all results remain very similar for 
small positive loads). At low temperatures the strongly 
fluctuating frictional force F x (t) (not shown) averages 
to a very small value, yielding a tiny friction coefficient 




FIG. 4: At T M a hard tip plowing through solid NaCl(100) is 
surrounded by a local liquid cloud that moves with it. Light 
blue spheres are tip atoms, green stand for Na + and yellow 
stand for CI - respectively. Note the fast healing furrow be- 
hind the tip. For clarity of illustration, the penetration here 
is 18 A. 
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FIG. 5: Frictional force of a grazing tip on NaCl(100). Note 
the very low friction, sharply increasing at Tm and above. 
Here the effective contact surface area was A -170 A 2 . 



fi = (F x )/F z ~ 0.007 at 300 K. Such nearly frictionless 
sliding of two hard incommensurate crystal lattices is not 
surprising, and resembles that recently reported for a ro- 
tated graphite flake on graphite 



As temperature 
grows, an initially moderate thermal growth of (F x ) up 
to 900 K gives way to an important nonlinear frictional 
surge near Tvi- This high temperature frictional peak 
makes a sharp contrast with the high temperature drop 
just seen in plowing friction on the same surface. Here 
the grazing frictional energy is dissipated mostly through 
generation of phonons in the solid substrate, plus some 
modest tip-induced drift of diffusive surface atoms taking 
place at the highest temperatures. Unlike low tempera- 
ture plowing, the solid surface is left, after dissipating 
away the small amount of frictional heat, very much in 
the same state it was before friction. This observation 
justifies an attempt at treating grazing friction within 
linear response theory, assuming a steady sliding state, 
and negligible influence of the tip on the surface. 
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Modeled as a hypothetical impurity-like particle, as- 
sumed free to move along (xy) on the hot surface, a tip 
will experience random forces - Brownian kicks. In the 
absence of an external force, it would undergo 2D diffu- 
sion; under a planar force /, it should drift with average 
velocity v proportional to /. Simulations at 900 K with 
increasing velocities between 15 and 50 m/s confirmed 
an approximately linear increase of friction force. Here 
the viscous behavior of grazing friction is due to ran- 
dom kicks from large thermal vibrations, rather than to 
viscoelastic properties of the hot surface, still perfectly 
solid and crystalline. Assuming a Langevin equation 
MR = —rjMR + f(t) for the in-plane Brownian mo- 
tion of the idealized (force-free) tip centered at R, the 
zero frequency self correlation of the random force f(t) 
is directly related to friction, in the simple form (lij 
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where it is moreover assumed, as usual that (/(<)) = 0. 
Within linear response, and for negligible tip influence 
over the surface, the random force acting on the tip can 
be taken to be proportional to the surface atom displace- 
ments Ui(t). Following Ying et al. [13, HH it is convenient 
to use 2D fc-space Ui(t) = J d 2 qu q (t)e~ lqRi and write 
the random tip force in the form f(t) = ~}2 q C q u q {t) 
where C q is a coupling function expressing the assumed 
proportionality between kicks onto the tip with individ- 
ual displacements Ui(t) and the random force f(t) felt 
by the tip. Inserting this expression into the equation for 
ij, we obtain 
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q \ 2 S{q,uj 



0) 



(2) 



where S{q, ui = 0) = J Q °° dtu^ q (t)u q (0) is the two dimen- 
sional dynamical structure factor of the tip-free surface 
substrate, proportional to the substrate atom displace- 
ment auto-correlation functions. In this manner, the 
grazing friction rj is reduced to a property of the unper- 
turbed surface. In the simple present form, it could also 
be arrived at directly by observing that the dissipation 
rate takes in linear response theory a Fermi Golden rule 
form, \C q \ 2 measuring the perturbation strength, and the 
dissipative part of the system's linear response function 
giving the density of final states. Considering that the 
frequencies that are relevant for dissipation (long wave- 
length phonons and diffusional processes) are extremely 
small compared to e.g. the solid substrate Debye fre- 
quency, one can set id = recovering (disregarding fac- 
tors) eq. ^ above. 

We can now attempt a comparison of the temperature 
dependence predicted by linear response theory eq. ([2|) 
of grazing friction based on the dynamical properties of 
the tip-free surface with the actual simulated friction. 



Here we deliberately oversimplify by assuming a T inde- 
pendent coupling C q - in reality, hard core interatomic 
potentials will generally imply a growing value with tem- 
perature, since vertical displacements, which eq. J2]) dis- 
regards, will be more important at higher temperature. 
In q space, C q will be nonzero for q < qo = 2w/L, where 
L is the lateral dimension of the tip (~ 13 A) or roughly 
go — 0.5 A -1 . Since go is small, we can simply assume 
Cq = Cq {qo — q), whence 
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To check the theory against frictional simulation, we 
extracted displacement auto-correlation functions of 
NaCl(100) from tip-free surface simulations. Fig.[6]shows 
S(lj) — J q<qo dqS(q,uj) as a function of temperature, of 
surface atoms only. Comparison with Fig. [5] is qualita- 
tively good, although as anticipated some increase of C q 
would be needed to make the agreement quantitative at 
high temperature. The sharp rise close to Tm and u — > 
is attributable to large and eventually (close to T ss ) catas- 
trophic softening of the very anharmonic surface lattice, 
still solid but extremely compliant in this regime. 

A well known analogue of this phenomenon is actu- 
ally realized in the physics of type II superconductors 
close to H C 2- Here the flux lattice - whose frictional de- 
pinning from the ion lattice and impurities determines 
the critical current - turns soft just before eventually 
disappearing j2^|. This soft state is so compliant that 
the pinning force, and the critical current with it, de- 
velops a last sharp peak before dropping to zero at i? C 2, 
a phenomenon that has been described by very similar 



formulas to those of sliding friction [2CJ, |21| . It is thus 
reasonable to propose that the high temperature friction 
increase of a hard grazing slider we just described is the 
direct frictional analogue of this peak effect in type II 
superconductors. 

In summary, the sliding of hard tips onto hard non- 
melting surfaces (here modeled by NaCl(100)) reveals 
new phenomena in high temperature nanofriction. As 
the melting point is approached, friction may either drop 
as in skating, or rise as in flux lattice depinning. We 
believe that these phenomena should be more general 
than the specific NaCl surface context in which they 
emerged. In particular, most metals possess at least 
one close packed non-melting surface, such as Pb, Al or 
Au(lll) Q, that should behave quite similarly to our 
results on NaCl(100). High temperature experiments on 
NaCl itself should be feasible. Because of the high vapor 
pressure (0.34 mmHg at Tm), and of a correspondingly 
high rate of evaporation (we estimate 1.7 x 10 ~ 5 s _1 , 
for terraces 50 A wide), the step flow velocity at sublima- 
tion conditions at Tm will be about 10 ~ 5 m/s. So long as 
the tip velocity is larger than this value, our description 
of nanofriction is fully applicable. Moreover, irrespective 
of sublimation the terraces between the steps will still 
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FIG. 6: Dynamical surface atom structure factor of NaCl 
S(uj — » 0) as a function of temperature for surface atoms 
only. Inset, S(ui) for temperature from 300 K to 1150 K. 



be well represented by the flat, dry, stable solid surface 
described in our simulations. Thus the effects described 
will be readily observed, for example once the sliding tip 
velocity v is faster than the step flow velocity. 
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